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Using the framework of soft-collinear effective theory (SCET), we factorize and calculate e+e^ 
angularity distributions, including perturbative resummation and the incorporation of a universal 
model for the nonperturbative soft function. Angularities T„ are a class of event shapes varying in 
their sensitivity to the substructure of jets in the final state, controlled by a continuous parameter 
fl < 2. We calculate the jet and soft functions in factorized angularity distributions for all a < 1 
to first order in the strong coupling a, and resum large logarithms to next-to-leading logarithmic 
(NLL) accuracy. We employ a universal model for the nonperturbative soft function with a gap 
parameter which cancels the renormalon ambiguity in the partonic soft function. 
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1. Introduction 

Event shapes yield simple information about the geometry of a hadronic final state in e^e~ 
annihilation and can be used to probe the strong interactions at various energy scales [1]. Two-jet 
event shapes e are designed so that they take a numerical value, usually between and 1, so that 
one of the kinematic endpoints (usually e = 0) corresponds to events with two perfectly-coUimated 
back-to-back jets in the final state. Event shape distributions depend on the hard-scattering cross- 
section at the large center-of-mass energy Q, on the perturbative branching and showering of the 
hard partons into jets at intermediate scales, and on the soft color exchange between jets and 
hadronization at a soft scale Aqcd- Event shapes are thus are particularly useful probes of both 
perturbative and nonperturbative effects in QCD, allowing, for instance, extraction of the strong 
coupling as and nonperturbative shape function parameters (see talks by I. Stewart and V. Mateu). 

In this talk, I overview the calculation, performed in [2] , of a particular class of event shape 
distributions, the angularities, using the framework of soft-colUnear effective theory (SCET) (see 
review talk by S. Fleming). We first establish a factorization theorem and calculate the hard, jet, and 
soft functions for angularity distributions to next-to-leading order (NLO) in a^. Then we solve the 
renormalization group equations for these functions and thereby resum large logarithms to next- 
to-leading logarithmic (NLL) accuracy in the angularity distributions near the two-jet kinematic 
endpoint. Finally we account for the effects of hadronization by introducing a universal model for 
the nonperturbative soft function, and present our final predictions for angularity distributions. 

2. Event Shapes and Angularities 

The most familiar event shape is the thrust, T = ^inaxtl^,gx |t"P(i' where Q is the e^e~ 
center-of-mass energy, and t, the thrust axis, is the unit three-vector which maximizes the sum 
of projections of final-state particles' three-momenta onto this axis. Once the thrust axis is deter- 
mined, many other event shapes can be defined, such as the jet broadening, B = ^ Jliex |t x P;!- A 
generalization of thrust and jet broadening is the class of angularities [3], 

^"(^) = ^Y^EiSm^'diil-cosdiy-^' = i £ |pf |e-l'^'l(i-') , (2.1) 
^ iex ^ iex 

where in the first form, £",■ is the energy of final-state particle / and 0,- is its angle with respect 
to the thrust axis. In the second form, pj is the ith particle's transverse momentum and T], its 
rapidity with respect to t. The parameter a can be any real number, — oo < a < 2 for Xa to be 
an infrared-safe observable. Two special cases are a = and a = I, which correspond to the 
thrust and jet broadening, Xq = I — T and Ti = B. It is known that the form of the factorization 
theorem which holds for the thrust distribution breaks down for the broadening distribution. By 
varying a between and 1, we can study how this factorization breaks down as a approaches 1 
[4]. More generally, being able to vary a provides a wealth of information about the final state 
that is otherwise obscured in looking at a single event shape in isolation. Roughly speaking, the 
angularity distributions for larger a are. dominated by jets which are very narrow, while those for 
smaller a —>■ —oo are dominated by jets which are wider. 
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3. Energy Flow and Factorization in SCET 

The angularity distribution for events e^e^ ^ X at center-of-mass energy Q is given by [5] 




d'x e'^- £ 4, (0| f{x)5{Za - |0) , (3.1) 

i=VA 



where q = (2,0), jyj^ = qTyj^q (with Ty = and = J^Js), and lJ^v are the leptonic tensors 
corresponding to these vector and axial currents (see [5]). We have made use of an operator 
which acts on states X and returns the eigenvalue Ta(X). This operator can be constructed from a 
momentum flow operator S'T{ri), which in turn is constructed from the energy-momentum tensor 
[5, 6, 7], 

S-ji^n) = =— / lim / dt niToi{t,Rn) , (3.2) 

cosh- Tj Jo Jo 

where n is a unit vector in the direction given by rapidity tj and azimuthal angle ^. /? is the radius 
of a sphere centered on the e^e^ collision. (S't{ti) acts on states X according to (^r(T]) \X) = 
Hiex |pf I ^{'H ~ 'Hi) \^)' where the transverse momentum pj and the rapidity T},- are measured with 
respect to the thrust axis.' In terms of (^^(tj), the operator can be constructed, 

1 poo 

Ta = - JT]e-l^l(i-«)^r(T]). (3.3) 

To factorize and evaluate the distribution Eq. (3.1) in SCET, we match the QCD current 
and the operator Tq onto operators in SCET. The current matches onto two-jet operators in SCET 
containing coUinear fields in two back-to-back directions [8, 9], 

= L L e'^''"-P---^-'C,^^^^^^^ , (3.4) 

where = T^jT^ = 7^/5, and the sums are over the direction n of the light-cone vectors n,n = 
(l,±n), and the label momenta p„ fj [10, 11]. The jet fields Xn.n = ^nn^n.n are built up from 
coUinear quark fields i^,,./^ and collinear Wilson lines W„ ,i, which have been decoupled from soft 
gluons through the BPS field redefinition with soft Wilson lines Yn,Yn [12]. The Wilson lines are 
defined in [10, 1 1]. The matching coefficient C„,i is calculable in perturbation theory. 

To match onto SCET operators, we simply replace the energy-momentum tensor To, in 
Eq. (3.2) with the energy-momentum tensor in SCET. Since the Lagrangian of SCET (after the 
BPS field redefinition) splits into separate purely collinear and purely soft parts, so does the energy- 
momentum tensor, and, thus, also the event shape operator. That is, f„ = f^' + f" + f*, where Tq'"'* 
are constructed as above, but using the energy-momentum tensor of only the «,«-collinear or soft 
Lagrangian of SCET. We can then factorize the angularity distributions Eq. (3.1) in SCET, 

~= H{Q-^) f dx: dxl dr^dixa - T« - - <)7;;(T-Myf (Tf , (3-5) 

ao dT„ J 

'The operator thus imphcitly depends on the thrust axis t of the final state. An operator which returns this axis 
was constructed in [5] and can be used here. In SCET, however, no such operator is required, as the thrust axis can be 
identified (up to power corrections) with the collinear jet direction n. This approximation breaks down for a > 1, which 
is one manifestation of the breakdown of factorization for these values of a [3, 5], 
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where Oq is the total e+e^ qq Born cross-section, the hard coefficient H is the squared ampli- 
tude of the two-jet matching coefficients, //( 2; /x) = \Cnn{Qn/2,—Qn/2;}x)\^, and the jet and soft 
functions Ja",Sa are given by matrix elements of coUinear and soft operators, 

d'x e''-Tr (0| Z«,eWa5(T;; - T^^.Q^p |0) (3.6) 

d'x e'^-Tr (0| Xf.,-Q{^)pS{z! - f:);t,,_e(0)« |0) (3.7) 

5„«;M) = ^Tr(0|Fl(0)7;(0)5«-f:)7„(0)F,(0)|0), (3.8) 
where the traces are over colors. 
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4. Renormalization Group Evolution and Resummation 

A fixed-order calculation of the angularity distributions will be divergent in the endpoint region 
Tq ^ where the event becomes more and more two-jet-like. The divergent terms ~ a" (In"' 'Ca)/Ta 
must be resummed to all orders in to yield a reliable, finite prediction. This resummation can be 
achieved by solving renormalization group equations for the hard, jet, and soft functions defined 
in the factorized angularity distributions above, and running each function from the scale at which 
large logarithms in it are minimized to the common factorization scale }X [13]. 

We begin by calculating the hard, jet, and soft functions to fixed order in perturbation theory 
using the Feynman rules of SCET [10, 12]. The result for the renormalized hard function is [9, 14] 



1 



2k 



IK^ 



+ ln^^ + 31n 
Q 



the renormalized jet function is [2] 



j:{x:-^i)=8{x:,)\\+ 
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71^ 3 — 5a + 



1-^V 4 12 \-a 
and the renormalized soft function is [2] 
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(4.1) 



(4.2) 



(4.3) 



(4.4) 



The superscript "PT" denotes that this is only the partonic soft function, calculable in perturbation 
theory. Later, we will convolute it with a model for the nonperturbative soft function. The pertur- 
bative calculations leading to these results fail to remain infrared-safe for a > 1 , manifesting the 
breakdown of factorization for these values of a [4] . 
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The hard function obeys the renormaUzation group (RG) equation 

Ai^//(!2;At) = 7//(MMe;At), 



where the anomalous dimension yu takes the form 



(4.5) 



(4.6) 



To first order in a^, Yii[(^s] = —10isCpl% and ///[oti] = — 3avCf/7r. Meanwhile, the jet and soft 
functions obey the slightly more complicated RG equation, 

d 



\i—F{x;ii)= Jjr'YF{r-T'-^)FiT';^), 



djJ. 

where F = J,S. The anomalous dimensions take the form 
7f(t-t';^) = -r/r[a,] 



1- 


"0(t-t')" 




z-z' 



In ^ 5(t - t') + 7f [a.v] 5(t - t') 



(4.7) 



(4.8) 



where jj = 2 — a and j's = 1, and to first order in a^, 
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The solution to the RG equation Eq. (4.5) for the hard function is 

and to the RG equation Eq. (4.7) for the jet and soft functions, 

F{x;n)= J dx'UF[T-T';pi,piQ)F[T';pio), 
where the evolution kernel Up is given by [15] 

UFir-'c';n,HQ) 
InEqs. (4.10) and (4.12), 
«f(At,Aio) 



rs[as]=0. (4.9) 



(4.10) 



r(-«f) 

da ^ , , 
jF JaAtio) j8[a] 

ccs{^l) da , , ^ 
7F[a]+2 





^ JfOf 


' 0(t-t') " 


(1 







da 



aApo) ^[«] 



rf[a] 



da' 



«,{/^o)i3[a'] 



(4.11) 
(4.12) 

(4.13) 
(4.14) 



/a,{W)) ^[«] 

where /3[a] is the beta function of QCD, and now F = H,J,S, with jn = 1- 

The solutions Eqs. (4.10) and (4. 1 1) for the hard, jet, and soft functions can be used to express 
these functions at the factorization scale /i in Eq. (3.5) in terms of their values at (arbitrary) scales 
l^H,j,s- These scales can be chosen to minimize the logarithms in the fixed-order hard, jet, and 
soft functions, and then the RG running of each function to the scale pL resums logarithms of 
jJL I \Ihj,s- After convoluting the functions to obtain the full distribution, we find that the initial 
scales should be chosen near \Lii ^ Q, \lj ^ Q'hi^^^ "\ and Hs ~ Q^a- The running to jJ. (typically 
chosen neai^ pi ^ Q) achieves resummation of the log"'Ta/TQ divergent terms in the fixed-order 
angularity distributions. In the following results, we used the anomalous dimension Thjs['^s\ to 
two loops and Yhj.s[(^s\ to one loop to achieve NLL accuracy in the resummation of these terms. 
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5. Model of Nonperturbative Effects 

A reliable prediction of the angularity distributions requires consistent incorporation of the 
above perturbative calculations and a universal model for the nonperturbative soft function. Such a 
model which interpolates consistently (i.e. with well-defined RG evolution) between perturbative 
and nonperturbative regimes is the convolution [16, 17], 

Q 

where is a model function based on that introduced for thrust and other event shape distribu- 
tions in [18], which we have generalized here to arbitrary [2]. is a "gap parameter" [16] that 
accounts for the minimum value of due to hadronization in each hemisphere. Specifically, 

2" 



Sa{Za,pi) = J dr'^sYir, - U " ^) , (5-1) 




/rP(T„) = 0(T«)^(A,5)|-(^) if, ( 1,1+A,^(^) )exp 



B+l [QTa 



(5.2) 

where A^, = A/(l — a) and A^ = A/(l — t?) are related to the a = (thi^ust) model parameters. We 
choose typical values for the model parameters (cf. [18]), A = 2.5, B = —0.4, and A = 0.55 GeV. 
The constant y^{A,B) normalizes the total integral of /^"^ to 1. Eq. (5.2) is the simplest general- 
ization of the thrust model function that obeys the universal scaling of the first moment of the soft 
function, {'Ca)s ~ 1/(1 — a), proposed in [19] and proven to all orders in a, in [7]. 

The gap parameter A^, and partonic soft function S^^ contain l/Q renormalon ambiguities, 
which can be removed by shifting the gap by a perturbatively-calculable amount, Aq = Aa{lJ.) + 
5a (/I ) where A^, is renormalon-free, and we choose the perturbative shift 5a to be 

2 /jT„e-eW(«er.)5PT(^^,.^) " 4^ ln^+£r^a,j, (5.3) 

where different choices of the free parameter R determine different renormalon subtraction schemes 
[20, 21]. Expanding the soft function Eq. (5.1) in powers of a^, S^^ combines with 5a -dependent 
terms to become renormalon-free, as is the new gap parameter A^. These schemes provide a 
consistent evolution equation for A^ in both n and R (see talk by I. Scimemi). Below, we choose 
the renormalon-free gap to be Aa{l^ = 1 GeV) = 100 MeV and R = 100 MeV. 



6. Final Predictions for Angularity Distributions 

In Fig. 1 we plot angularity distributions for —2 < a < 1/2 for 2 = 100 GeV, including 
the above ^(a^) fixed-order calculations of the hard, jet, and soft functions, the NLL resumma- 
tion in the Tq ^ endpoint region to all orders in as, and the soft model function of Eq. (5.2). 
We also matched the resummed distribution, which is accurate in the endpoint region, to fixed- 
order QCD, which is more accurate in the tail region, using the procedure described in [2]. We 
choose the hard scale to be Hh = Q, and the jet and soft scales to interpolate between a mini- 
mum fixed value 5 = /x^^" > Aqcd for small Za and jXj = Qr}/^^ and /I5 = Qta for larger 
Tfl, allowing us to avoid spurious Landau poles [2]. In the plots, we chose /i^'" = 1 GeV and 
^min ^ g(i-.,)/(2-«)(^min)i/(2-a) pjg 2 we plot the a = -1,0,1/2 distributions varying the 
scales Hj s together (coiTclated) by factors of 2. Dependence on /i and Hh is much weaker. 
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Ta 

Figure 1: Angularity distributions for —2 < a < j at Q = 100 GeV, with ^{(Xs) hard, jet, and soft functions, 
NLL resummation, and gapped model soft function. 

1 da 
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Figure 2: Correlated variation of jet and soft scales pLj and jis together by factors of 1 /2 and 2 in angularity 
distributions at 2 = 100 GeV for a= — 1, a = 0, and a = 

7. Conclusion 

We have predicted angularity distributions resummed in perturbation theory to NLL accuracy, 
including for the first time the jet and soft functions in the factorization theorem to NLO, and a 
universal model for the nonperturbative soft function. We used the framework of SCET to perform 
the factorization, resummation, and incorporation of the nonperturbative model in a unified way. 
Comparison to data from LEP or a future linear collider will test the robustness of the model 
we employed for the nonperturbative soft function. Extension of the notion of event shapes to 
individual "jet shapes" can also allow the probing of jet substructure in hadron collisions [22, 23]. 

This work was supported in part by the U.S. Department of Energy under Contract DE-AC02- 
05CH11231 and the National Science Foundation under grant PHY-0457315. 
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